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A floating Wigner crystal differs from the standard one by a spatial averaging over positions of 
the Wigner-crystal lattice. It has the same internal structure as the fixed crystal, but contrary to it, 
takes into account rotational and/or translational symmetry of the underlying jellium background. 
We study properties of a floating Wigner molecule in few-electron spin-polarized quantum dots, and 
show that the floating solid has the lower energy than the standard Wigner crystal with fixed lattice 
points. We also argue that internal rotational symmetry of individual dots can be broken in arrays 
of quantum dots, due to degenerate ground states and inter-dot Coulomb coupling. 

PACS numbers: 73.20.Qt, 73. 21. La, 73.22.Gk 

At sufficiently low densities a system of interacting electrons on a uniform positive background forms a correlated 
Wigner-crystal state. The idea of an electron crystal was originally proposed by Wignerjd and then extensively 
studiccLiii three-tnia and two-tdimensionaltjCJ (2D) electron systems (experimental realizations are electrons on liquid 
HeliumEa and semiconductorc^l surfaces). Recently, an interest aroused to the formation and properties of Wiener 
molecules in semiconductor quantum dots - zero-dimensional systems with a finite number N of 2D electrons .EZIlj 

A standard Wigner-crystal ground-state trial wave function ^wc{ r \i r 2, ■ ■ ■ ,^n) (for a fully-spin-polarized elec- 
tron system) has the form of a Slater determinant constructed from single-particle orbitals ip( r ) centered at points, 
corresponding to equilibrium positions of classical particles, 

*wc({ri}, {R,}) = det |^(r 4 - Rj-)| . (1) 

In a circular parabolic quantum dot with a small numbe^JV (< 8) of 2D electrons the vectors Rj form symmetric 
polygonal configurations with N c electrons in the center f3t3 

Rj = i?(cos^,sin^), 6j = 2n(j - 1)/N, j = l,...,N, N < 5, 
Ri=0, Rj = E(cos0j,sin0j), % = 2tt(j - 2)/(JV - 1), j = 2,...,N, N = 6,7,8, (2) 

N c = at N <_5u_and N c — lat6<N<8 (for quantum dots with more particles see exajpcles of classical 
configurations inEZlEa). In a macroscopic 2D electron system vectors Rj form a triangular lattice,li3liZl 

{Rj} ~a(h + h/2,V3l2/2), n s a 2 = 2/V3, (3) 

with the period a corresponding to the average area density of 2D electrons n s (li and I2 are integer). 

The equilibrium classical configurations (|^) are usually considered to be fixed with respect to some fixed external 
reference frame. This leads, in fact, to a violation of the rotational (in circular quantum dots) and translational (in 
infinite 2D systems) symmetry of the ground state, and results in an angular- and position-dependent ground-state 
electron density. The trial many-body wave function ([!]) in the dots (in an infinite 2D electron system) is not the 
eigenfunction of the total angular momentum L to t (of the total momentum Ptot), although the Hamiltonian has the 
rotational (translational) symmetry and commutes with the corresponding operators (here L to t = where z-axis 
is perpendicular to the plane of the 2D system). 

The symmetries can be restored in (Q) by averaging over the corresponding positions of the v gGtO |rs Rj. Such a 
floating solid was discussed in the literature in connection with other many-body problems, see e.g.rlr 8 !. In this paper 
we discuss some general properties of a floating Wigner molecule in a circular parabolic quantum dot, and of a floating 
Wigner crystal in a macroscopic 2D electron system. We show that {R^ j-averaged Wigner-crystal wave functions can 
be constructed as eigenfunctions of the total angular/total momentum operators, and that these eigenstates possess 
the symmetry of the Hamiltonian. We perform specific calculations for the energy of few-electron quantum dots, in 
the fixed and floating Wigner-crystal states, and show that the {R 3 ;}-averaging does lead to an essential reduction 
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of the variational ground-state energy. For two-, three- and four-electron parabolic quantum dots we compare our 
variational results with exact results available in the literature. We also show that the internal rotational symmetry 
of the wave functions in circular dots can be broken if the ground state is degenerate with respect to the total angular 
momentum L: a superposition of eigenfunctions with different L breaks the symmetry without increasing the energy of 
the state. We argue that this leads to phase transitions in quantum-dot arrays with spontaneous symmetry breaking. 

Consider a circular TV-electron parabolic quantum dot. Let R" = Rj(cos(6j + a),s'm(9j + a)) be a set of vectors 
rotated around the origin by an angle a with respect to (Q). Consider a function 

= f ^ LQ ^c({rJ, {R"}) 5 (4) 
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obtained from ([!]) by averaging over a with the weight function e lLa , and labeled by an arbitrary integer L. The 
functions (|J) have the same internal structure as (l]), so that electron-electron correlations are taken into account at 
the same level of accuracy, and satisfy the eigenvalue equation 

L tot ^ L ({Ti}) = hL^ L ({n}), (5) 

with L being the total angular momentum quantum number. The fixed Wigner-crystal wave function (|l]) can be 
presented in the form of an infinite sum of eigenfunctions ^l with all possible values of L, 

oo 

9 wc ({Ti},{*i})= E (6) 

L — — oo 

These results obviously remain valid if the trial wave functions are multiplied by additional Jastrow factors (see 
e.g.c3c3) which depend only on relative coordinates of electrons. 

For a quantum dot with N < 8 electrons we calculate expectation values E = (^\H\^!) / (^f\^f) of the Hamiltonian 
H in the states (f|) and (Q), and directly compare the energies El and Ewe of the floating and fixed Wigncr- molecule 
states. The Hamiltonian of the dots has the form 



2m* 2 ) 2 i? , 

where m* is the effective mass of 2D electrons, k is the dielectric constant of the host semiconductor, and loq charac- 
terizes the confinement. For the single-particle orbitals ip(r) we take gaussians 

^ r) -^ eXP ("S)' (8) 

with the width I being either equal to the oscillator length Iq = \f%[m~*ujQ, or considered as a variational parameter. 
For the vectors Rj we assume the classical equilibrium configurations (^|). 

Due to the symmetry of configurations (||) not all values of L are allowed in Eq. (^). For instance, if N = 2 and the 
system is fully spin-polarized (the total spin is S = N/2 = 1) the number L can take only odd values; otherwise the 
integral (Q) vanishes. For the ground state it is reasonable to assume the smallest possible value of the total angular 
momentum, so that we get L — ±1. As known from exact-diagonalization calculations,!^ at S — 1 the ground state 
of a two-electron quantum dot has indeed the total angular momentum \L\ = 1. Similarly, if N = 3, allowed values of 
L are 0,±3,±6, . . ., and in the ground state of the-£ully polarized (S = 3/2) three-electron dot we have L = 0. This 
also agrees with exact-diagonalization calculationsL3. In general, for iV-electron dots with N from 2 to 8, the total 
spin S = N/2, and classical configurations (0), we get the following ground-state total angular momenta 

L=(N- N c )/2, if (N - N c ) = even, 

L=Q, if (AT — N c ) = odd, (9) 

where N c = 0, if N < 5, and N c = 1, if 6 < N < 8. 

The energies El (with L from Eq. @) and Ewe have been calculated with the help of general formulas obtained 
in Ref.c2l. Figure [j] shows the difference (Ewe — E exact ), (El — E cxact ) between the,-£iiergies of the fixed and 
floating variational Wigner-crystal states and exact-diagonalization results from RefsOOllj, for N = 2, 3, and 4 
spin-polarized electrons, S = N/2. Figure || exhibits the difference (Ewe ~ El) between the energies of the fixed 
and floating Wigner-molecule states, for parabolic quantum dots with N from 2 to 8 electrons. One sees that in all 
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the cases the Rj-averaging leads to an essential improvement (by about 0.2-0.3 Tiuiq) of the variational ground state 
energy, and that in general the energy gain increases with the number of particles in the dots. Another (quite natural) 
finding is that for "J l states with zero total angular momentum (four upper thick curves in Figure ||) the energy gain 
is substantially larger than for the states with L ^ (three lower thin curves). The difference between the variational 
and exact ground-state energies is reduced by a factor of w 2 in the two-, three-, and four-electron quantum dots, 
Figure [l], due to the floating version of the Wigner-molecule trial wave function. 

The states are eigenstates of the total angular momentum Ltot, Eq. (||), and the density of electrons n^(r) in 
these states does not depend on the angular coordinate, til{t,0) = til{t)- However, if the ground state is degenerate, 
as it is the case, for instance, at L ^ (N = 2, 4, and 7), one can construct solutions 

*mixed = C L V L + C*_ L *_ L , (10) 

which have the same (ground-state) energy E m ixed = E±l, but an angular-dependent density n m i xec i(r, 9). A possible 
mixed state is (^l + ^-i)/v2 with (L) — 0. Such broken-symmetry ground states can thus exist in a circularly 
symmetric quantum dot, depending on initial conditions. In a single dot the energies of symmetric and broken- 
symmetric solutions are the same, but in the presence of a small perturbation the symmetry-broken state may have 
the lower energy (an evident example is a circular dot with "disorder" , e.g. with an asymmetrically located impurity). 
Consider for instance an array of circular dots with a weak inter-dot Coulomb coupling, and assume that each dot 
contains two spin-polarized electrons (no tunneling is assumed between the dots). Then the structure schematically 
shown in Figure^ will obviously have the lower energy than the state with angular independent density in each dot. A 
weak inter-dot Coulomb interaction in arrays of dots can thus lead to phase transitions with a spontaneously broken 
phase. 

Results presented in formulas (Q) - (||) can be easily generalized to the case of other electron systems on the jellium 
background. For example, for an infinite 2D electron system (with periodic boundary conditions) the floating Wigner 
crystal wave function can be written as 

,iKa 



*K({r<» = J dae^tfwc ({ri}, {R,- + a}), (11) 

where the integral is taken over the (large) area of a sample. The functions (^Tj) are the eigenfunctions of P tot , 

Pfot*K({r<» = ^K* K ({rJ), (12) 

with hK. being the total momentum quantum number. The fixed Wigner crystal wave function (Q) is expanded in 
terms of as 

*wc{{Ti}, {R;}) = j-^ J dK* K ({ri». (13) 

The floating Wigner crystal wave functions ( |ll| ) give the position-independent density of 2D electrons. However, if 
the ground state is degenerate with respect to K, a superposition of such degenerate ground states can again lead to 
a modulated density. 

To summarize, we have studied some general properties of a floating Wigner crystal in parabolic quantum dots and 
other electron systems. In few-electron dots the rotationally-invariant floating solid was proved to have a substantially 
lower energy than the fixed Wigner molecule given by the trial wave function of Eq. (|l|) . In arrays of quantum dots 
the inter-dot interaction may lead to a phase transition to the ground state with broken rotational symmetry in 
individual dots. 

The work was supported by the Deutsche Forschungsgemeinschaft through the SFB 484. We thank Karl-Heinz 
Hock for useful discussions. 
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(curves labeled as "fixed") and El — E exac t (curves labeled as "floating") for 
2, (b) N = 3 and (c) TV = 4 spin-polarized electrons, as a function of the Coulomb 



FIG. 1. Energy difference Ewe - 
parabolic quantum dots with (a) TV 
interaction parameter lo/as = v/ e 2 /aBffjJo (energy unit is TtoJn, as is the effective Bohr radius). For both trial wave functions 
the curves without and with optimization over the variational parameter I are shown tihe curves labeled as, "/ = in" and "Z 
varied" respectively). Exact-diagonalization results are taken from Ref.L3 (N — 2), Ref.EJ (N = 3), and RefO (N = 4). 
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FIG. 2. Energy difference Ewe — El, in units hu>o, as a function of the interaction parameter Io/cib, for parabolic quantum 
dots with N = 2 to 8. For both trial wave functions optimization over the variational parameter I have been performed. 




FIG. 3. Possible phase transition to a broken-symmetry state in an array of Coulomb-coupled quantum dots. It is assumed 
that there are two electrons in each dot (thin circles show the maxima of the electron density), and that electron spins are 
polarized. 
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